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Abstract — In view of the discovery of a new boson by the ATLAS and CMS Collaborations at the LHC, 
we present an update of the global Standard Model (SM) fit to electroweak precision data. Assuming the 
new particle to be the SM Higgs boson, all fundamental parameters of the SM are known allowing, for the 
first time, to overconstrain the SM at the electroweak scale and assert its validity. Including the effects of 
radiative corrections and the experimental and theoretical uncertainties, the global fit exhibits a p-value 
of 0.07. The mass measurements by ATLAS and CMS agree within 1.3er with the indirect determination 
M H = 94 ±22 GeV. Within the SM the W boson mass and the effective weak mixing angle can be accurately 
predicted to be M w = 80.359±0.011GeV and sin 2 ^ ff = 0.23150±0.00010 from the global fit. These results 
are compatible with, and exceed in precision, the direct measurements. For the indirect determination of 
the top quark mass we find m t — 175.8 ± 2 ' 4 GeV, in agreement with the kinematic and cross-section based 
measurements . 
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1 Introduction 

The discovery by the ATLAS [1] and CMS [2] experiments at the LHC of a new particle with 
mass ~126 GeV and with properties compatible with those of the Standard Model (SM) Higgs 
boson concludes decades of intense experimental and theoretical work to uncover the mechanism 
of electroweak symmetry breaking and mass generation. If forthcoming data confirm that the 
new particle is the SM Higgs boson, this discovery exhibits another - possibly the greatest ever - 
triumph of the SM, as not only the SM predicts the Higgs couplings to the SM fermions and bosons, 
but it also constrains the Higgs boson to be light compared to its unitarity bound of roughly a TeV. 
This indirect information on the Higgs mass was extracted from Higgs loops affecting the values 
of Z boson asymmetry observables and the W mass. Global fits to precisely measured electroweak 
data derived 95% confidence level (CL) upper limits on the Higgs mass of around 160 GeV [3-6]. 

In this letter we interpret the new particle as the SM Higgs boson and present the consequences 
on the global electroweak fit. A detailed description of the experimental data, the theoretical 
calculations, and the statistical framework used in the analysis is provided in past publications [6, 
7]. Here, we only briefly recall the most relevant aspects of the analysis and highlight recent 
changes. The main goal of this letter is to quantify the compatibility of the mass of the discovered 
boson with the electroweak precision data and its impact on the indirect determination of the W 
boson mass, the effective weak mixing angle, and the top quark mass. The implications of the 
discovery on the SM with three and four fermion generations were also studied in [8]. 

2 Experimental data and theoretical predictions 

The experimental inputs used in the fit include the electroweak precision data measured at the 
Z pole and their correlations [9], the latest world average values for the mass of the W boson, 
M w = 80.385±0.015GeV [10], and its width, T w = 2.085±0.042GeV [11], the latest average of the 
direct top mass measurements from the Tevatron experiments, mt = 173.18 ± 0.94 GeV [12], 1 and 
the hadronic contribution to the running of the electromagnetic coupling strength, Aa{^ d (M|) = 
(2757 ± 10) • 10 -5 [19]. For the Higgs boson mass we use the measurements from ATLAS, Mjj = 
126.0 ± 0.4 ± 0.4 GeV [1], and CMS, M H = 125.3 ± 0.4 ± 0.5 GeV [2], where the first uncertainties 
are statistical and the second systematic. A detailed list of all the observables, their values and 
uncertainties as used in the fit, is given in the first two columns of Table 1. 

A proper average of the Higgs mass measurements requires a detailed experimental study of the 
systematic correlations. Owing to the weak (logarithmic) dependence of the electroweak fit on 
the Higgs mass, we find that the fit results are insensitive to the difference between a straight 
(uncorrelated) weighted average of the ATLAS and CMS measurements {Mh = 125.7 ± 0.4 GeV), 

1 The theoretical uncertainties arising from nonperturbative colour-reconnection effects in the fragmentation pro- 
cess [13, 14], and from ambiguities in the top-mass definition [15, 16], affect the (kinematic) top mass measurement. 
Their quantitative estimate is difficult and may reach roughly 0.5 GeV each, where the systematic error due to shower 
effects could be larger [13]. To estimate the effect of a theoretical uncertainty of 0.5 GeV, inserted in the fit as a 
uniform likelihood function according to the _Rfit scheme [17, 18], we have repeated the indirect determination of 
some of the most relevant observables. We find in particular Mh = 90 IJi GeV, Mw = 80.359 ± 0.013 GeV, and 
sin 2 #f ff = 0.23148 ± 0.00010, which, compared to the standard results given in Eq. (1), (4) and (7), exhibit only a 
small deterioration in precision. 
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and their weighted average obtained assuming the systematic uncertainties to be fully correlated 
(Mji = 125.7 ± 0.5 GeV). In this paper the former combination is used. 2 

For the theoretical predictions, we use the calculations detailed in [6] and references therein. They 
feature among others the complete O(a^) calculation of the QCD Adler function [20, 21] and the 
full two-loop and leading beyond-two-loop prediction of the W mass and the effective weak mixing 
angle [22-24]. Two modifications apply here: first, an improved prediction of R® is invoked that 
includes the calculation of the complete fermionic electroweak two-loop (NNLO) corrections based 
on numerical Mellin-Barnes integrals [25]; second, the calculation of the vector and axial- vector 
couplings, g^ A and g(r, now entirely relies on accurate parametrisations [26-29]; the correction 
factors from a comparison with the Fortran ZFITTER package [30, 31] applied previously at very 
high Higgs masses [6] are no longer used. 

3 Results 

The fit to all input data from Table 1 converges with a global minimum value for the test statistics 
of Xmin = 21-8, obtained for 14 degrees of freedom. Using a pseudo Monte Carlo (MC) simulation 
and the statistical method described in [7] we find the x m in distribution shown in Fig. 1. The 
resulting p- value for the SM to describe the data amounts to 0.07 (corresponding to 1.8cr). This 
result is consistent with the naive p-value Prob(21.8, 14) = 0.08. 

The inferior compatibility of the fit compared to earlier results [6] is not primarily caused by the 
insertion of the new Mh measurements in the fit, but is due to the usage of a more accurate R® 
calculation [25] that leads to a smaller SM prediction. 3 

The results of the complete fit for each fit parameter and observable are given in the fourth column 
of Table 1, together with their uncertainties estimated from their A% 2 = 1 profiles. Figure 2 (left) 
shows the pull values obtained from the difference between the result of the fit and the input data 
in units of the data uncertainty. No single pull value exceeds 3a. The known tension between the 
left-right asymmetry and Ap^ is reproduced. The new R® calculation [25] increases the discrepancy 
between the R® prediction and its measurement from —0.8a to — 2.4cr. 

The fifth column in Table 1 gives the results obtained without using the Mh measurements in 
the fit (i.e., Mg is a freely varying parameter). In this case, which represents the well known 
result of the standard electroweak fit prior to the Mh measurement, the fit converges with a global 
minimum of Xmin = 20.3 for 13 de grees of freedom (Prob(20.3, 13) = 0.09). We obtain 

M H = 94 ±! GeV . (!) 

consistent within 1.3c with the Mh measurements. The top left panel of Fig. 3 displays the 
corresponding A% 2 profile versus Mh (grey band) compared to the new Mh measurements of 
ATLAS and CMS (red/orange data points) and the A% 2 profile of the fit including the Mh 
measurement (blue curve). 

2 The main source of systematic uncertainty in the ATLAS and CMS mass measurements stems from the energy 
and momentum calibrations, which should be uncorrelated between the experiments. 
3 The quantity of 7?° has only little dependence on Mh [25]. 
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Parameter 


Input value 


Free 
in fit 


Fit result 
incl. Mh 


Fit result 
not incl. Mh 


Fit result incl. Mh 
but not exp. input in row 


l\lH [OeVp ' 


Izo. ( ± U.4 


yes 


1 0£ 7 I n 4 
tlD.t ± U.4 


94 -22 


QA + 25 
94 -22 


M w [GeV] 


80.385 ±0.015 


- 


80.367 ±0.007 


80.380 ±0.012 


80.359 ±0.011 


T w [GeV] 


2.085 ± 0.042 


- 


2.091 ±0.001 


2.092 ±0.001 


2.091 ± 0.001 


My [GeVl 


91.1875 ± 0.0021 


ves 


91.1878 ± 0.0021 


91.1874 ±0.0021 


91.1983 ± 0.0116 


Ty [GeVl 


2.4952 ± 0.0023 




2.4954 ± 0.0014 


2.4958 ± 0.0015 


2.4951 ± 0.0017 


"had L 11LJ J 


41.540 ± 0.037 




41.479 ± 0.014 


41.478 ± 0.014 


41.470 ± 0.015 


#° 


20.767 ± 0.025 




20.740 ± 0.017 


20.743 ± 0.018 


20.716 ± 0.026 


A o,i 

^FB 


0.0171 ±0.0010 


_ 


0.01627 ± 0.0002 


0.01637 ± 0.0002 


0.01624 ± 0.0002 


A e W 


0.1499 ±0.0018 


- 


n 1/17Q +0.0006 

u.i4/o_ 0008 


0.1477 ±0.0009 


0.1468 ±0.0005 (t) 


sin 2 ^ ff (Q FB ) 


0.2324 ±0.0012 


- 


n oqi as +0.00011 
U.zdl48_ 00007 


n ooi /I q +0.00010 
U.ZOI^O _ 00012 


0.23150 ± 0.00009 


A c 


0.670 ±0.027 


- 


o-668o±8:8gggS 


n «fioo +0.00042 

U.bb8z_ 00035 


0.6680 ± 0.00031 


-4 (J 


0.923 ±0.020 


- 


0.93464 +£888$ 


0.93468 ± 0.00008 


0.93463 ± 0.00006 


,0,c 
^FB 


0.0707 ±0.0035 


- 


073Q +0 0003 
u - u ' OJ -0.0005 


0.0740 ± 0.0005 


0.0738 ± 0.0004 


^FB 


0.0992 ± 0.0016 




0-1032 +°;°°°* 


0.1036 ± 0.0007 


0.1034 ± 0.0004 




u.i i zi m u.uuou 




U. 1 iZZO zc u.uuuuu 


u. i * zzo m u.uuuuo 


n 1 799'? + n nnnnfi 

U.I 1 ZZO m U.UUUUO 




0.21629 ± 0.00066 




0.21474 ± 0.00003 


0.21475 ± 0.00003 


0.21473 ± 0.00003 


rn c [GeV] 


i o7 +0.07 
1 - z ' -0.11 


yes 


i 07 +0.07 
-0.11 


i 07 +0.07 
1 - Z( -0.11 




m b [GeV] 


4 20 + 17 

^- zu -0.07 


yes 


4 20 +0 ' 17 

^■ zu -0.07 


4 20 + 17 

^■ zu -0.07 




m t [GeV] 


173.18 ±0.94 


yes 


173.52 ±0.88 


173.14 ±0.93 


175.8±2:I 


A^d^i) (Av) 


2757 ± 10 


yes 


2755 ± 11 


2757 ± 11 


2716 ±% 


a s (M|) 




yes 


0.1191 ±0.0028 


0.1192 ±0.0028 


0.1191 ±0.0028 


<5thA% [MeV] 


[-4,4]thco 


yes 


4 


4 




<5 th sin 2 ^ ff ( A > 


[-4.7,4.7] thco 


yes 


-1.4 


4.7 





(o) Average of ATLAS (Mh = 126.0 ± 0.4 (stat) ± 0.4 (sys)) and CMS (M H = 125.3 ± 0.4 (stat) ± 0.5 (sys)) 
measurements assuming no correlation of the systematic uncertainties (see discussion in Sect. 2). Average of 
LEP (Ae = 0.1465 ± 0.0033) and SLD (Ae = 0.1513 ± 0.0021) measurements, used as two measurements in the fit. 
(t) The fit w/o the LEP (SLD) measurement gives At = 0.1474 l°:Sool ( A e = °- 1467 -0.0004) • 
^'in units of 10 -5 . ' v 'Rescaled due to as dependency. 



Table 1 : Input values and fit results for the observables and parameters of the global clcctroweak fit. The 
first and second columns list respectively the observablcs/paramctcrs used in the fit, and their experimental 
values or phenomenological estimates (see text for references) . The subscript "theo" labels theoretical error 
ranges. The third column indicates whether a parameter is floating in the fit. The fourth column quotes 
the results of the complete fit including all experimental data. The fifth column gives the fit results for 
each parameter without using the Mjj measurement in the fit. In the last column the fit results are given 
without using the corresponding experimental or phenomenological estimate in the given row. 
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Figure 1 : Result of a pseudo MC analysis of the complete electroweak fit. Shown are distributions of the 
Xmin test statistics obtained from a pseudo MC simulation with varying (hatched histogram) and fixed 
theoretical uncertainty parameters (5th) in the fit (shaded/green histogram). The Xmin obtained in the 
complete fit to data is indicted by the arrow together with the p-values found for these two cases. Also 
shown is an idealised x 2 function assuming a Gaussian case with 14 degrees of freedom (solid black line) . 

Figure 2 (right) shows the determination of Mjj in fits in which among the four observables 
providing the strongest constraint on M#, namely A^(LEP), ^(SLD), ^4p B , and Mw, only the 
one indicated in a given row of the plot is included in the fit. For comparison also the indirect fit 
result (without the Mjj measurement) and the direct measurement are shown as vertical bands. 

The remaining plots in Fig. 3 show the Ax 2 profile curves versus mt (top right), M\y (bottom left), 
and sin 2 #^ ff (bottom right) obtained without using the corresponding experimental measurement in 
the fit (indirect determination, cf. last column of Table l). 4 For comparison also the corresponding 
profile curves excluding in addition the new Mjj measurements are shown (grey bands). The results 
from the direct measurements for each variable are also indicated by data points at Ax 2 = 1. 
The insertion of Mh substantially improves the precision of the fit predictions. 

4 For the indirect determination of sin 2 #f ff , shown as the blue band in Fig. 3 (bottom right), we exclude from the 
fit all experimental measurements with direct sensitivity to sin 2 #f ff , namely the measurements of Yz, o"h a d> Ri, ^fb> 
At, sin 2 ^ ff (Qpb), A c , Ab, A 1 ^^, -A™, B!i and R®. As a compensation of the missing value of R® we provide a value 
for a s (Mz). Since the fit results are independent of the exact a s value, we use our fit result 0.1191 ± 0.0028 in this 
case. 

5 We show the aforementioned result of the Tevatron combination of the direct top mass measurements [12], the 
top pole mass derived from the measured ti cross-section at the Tevatron (m t = 173.3 ± 2.8 GeV), assuming no new 
physics contributes to this cross section measurement [32], the direct top mass measurement of ATLAS determined 
in 1.04fb _1 of pp collisions at y/s = 7TeV (m t = 174.5 ± 2.4GeV) [33], the direct top mass measurement of CMS 
based on 5.0fb _1 of 7TeV data (m t = 173. 5± 1.1 GeV) [34], the aforementioned W mass world average [10] and the 
LEP/SLD average of the effective weak mixing angle (sm 2 9 e Eft = 0.23153 ± 0.00016) [9]. 
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Figure 2: Left: pull comparison of the fit results with the direct measurements in units of the experimental 
uncertainty. Right: determination of Mh excluding the direct Mh measurements and all the sensitive 
observables from the fit, except the one given. Note that the fit results shown are not independent. 



The fit indirectly determines the W mass (cf. Fig. 3 - bottom left, blue band) to be 

M w = 80.3593 ±0.0056 mt ±0.0026A/ z ±0.0018A Qhad (2) 
± 0.0017 as ± 0.0002 Mh ± 0.0040 thc o , (3) 

= 80.359 ± O.Olltot , (4) 

which exceeds the experimental world average in precision. The different uncertainty contribu- 
tions originate from the uncertainties in the input values of the fit as given in the second column 
in Table 1. The dominant uncertainty is due to the top quark mass. Due to the weak, logarith- 
mic dependence on Mh the contribution from the uncertainty on the Higgs mass is very small 
compared to the other sources of uncertainty. Note that in the i?fit scheme [17, 18] the treatment 
of the theoretical uncertainty as uniform likelihood corresponds a linear addition of theoretical 
and experimental uncertainties. Quadratic addition would give a total uncertainty in the M\y 
prediction of 0.008. 

The indirect determination of the effective weak mixing angle (cf. Fig. 3 - bottom right, blue 
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Figure 3: Ax 2 profiles as a function of the Higgs mass (top left), the top quark mass (top right), the W 
boson mass (bottom left) and the effective weak mixing angle (bottom right). The data points placed along 
Ax 2 = 1 represent direct measurements of the respective observable and their ±1<t uncertainties. The grey 
(blue) bands show the results when excluding (including) the new Mh measurements from (in) the fits. 
For the blue bands as a function of m t , Mw and sin 2 6r ff the direct measurements of the observable have 
been excluded from the fit in addition (indirect determination). The solid black curves in the lower plots 
represent the SM prediction for sin 2 6^ ff and M w derived from the minimal set of input measurements, as 
described in the text. In all figures the solid (dotted) lines illustrate the fit results including (ignoring) 
theoretical uncertainties in the fit. 



band) gives 

sin 2 (9ffT = 0.231496 ± 0.000030 mt ± 0.000015 A / Z ± 0.000035 A a had (5) 
± 0.000010 as ± 0.000002 Mh ± 0.000047 t heo , (6) 

= 0.23150 ± 0.00010 to t , (7) 

which is compatible and more precise than the average of the LEP/SLD measurements [9]. The 
total uncertainty is dominated by that from Aa%ad an d mt, while the contribution from the uncer- 
tainty in Mh is again very small. Adding quadratically theoretical and experimantal uncertainties 
would lead to a total uncertainty in the sin 2 i9^ ff prediction of 0.00007. 

Finally, the top quark mass, cf. Fig. 3 (top right, blue band), is indirectly determined to be 

m t = 175.8 ±2A GeV > ( 8 ) 
in agreement with the direct measurement and cross-section based determination (cf. Footnote 5). 



3 Results 



8 




Figure 4: Contours of 68% and 95% CL obtained from scans of fixed Myy and m t . The blue (grey) areas 
illustrate the ht results when including (excluding) the new Mh measurements. The direct measurements 
of M~w and m t are always excluded in the ht. The vertical and horizontal bands (green) indicate the la 
regions of the direct measurements. 

The measured value of Mh together with the fermion masses, the strong coupling strength cv s (M|) 
and the three parameters defining the electroweak sector and its radiative corrections (chosen 
here to be Mz, Gf and Aa^ d (M^)) form a minimal set of parameters allowing one, for the 
first time, to predict all the other SM parameters/observables. A fit using only this minimal 
set of input measurements 6 yields the SM predictions My/ = 80.360 ± 0.011 GeV and sin 2 ^ ff = 
0.23152 ± 0.00010. The A X 2 profile curves of these predictions are shown by the solid black lines 
in Fig. 3 (bottom left) and (bottom right). The agreement in central value and precision of these 
results with those from Eq. (4) and (7) (cf. blue bands in the plots) illustrates the marginal 
additional information provided by the other observables. 

Figure 4 displays CL contours of scans with fixed values of My/ and mt, where the direct measure- 
ments of My/ and mt were excluded from the fit. The contours show agreement between the direct 
measurements (green bands and data point), the fit results using all data except the My/, mt and 
Mh measurements (grey contour areas), and the fit results using all data except the experimental 
My/ and wi[ measurements (blue contour areas). The observed agreement again demonstrates the 
impressive consistency of the SM. 

Following the approach in [6] we extract from the electroweak fit the S, T, U parameters [35, 36] 
describing the difference between the oblique vacuum corrections as determined from the experi- 
mental data and the corrections expected in a reference SM (SM re f defined by fixing mt and Mh)- 
After the recent discovery, we change our definition of the reference SM for the S, T, U calculation 



6 For a s (M|) we use the result from Table 1. 
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Figure 5: Experimental constraints on the S and T parameters with respect to the SM reference (Mn.icf — 
126 GeV and m t , TC { = 173 GeV). Shown are the 68%, 95% and 99% CL allowed regions, where the third 
parameter U is left unconstrained (orange, left) or fixed to (blue, right). The prediction in the SM is given 
by the black (grey) area when including (excluding) the new Mh measurements. 



to M# jre f = 126 GeV and m t:Te f = 173 GeV. With these we find: 

S = 0.03 ±0.10, T = 0.05 ±0.12, U = 0.03 ± 0.10 , (9) 

with correlation coefficients of ±0.89 between S and T, and —0.54 (—0.83) between S and U (T 
and U). Fixing U = we obtain S\ u=0 = 0.05 ± 0.09 and T\ u=0 = 0.08 ± 0.07 with a correlation 
coefficient of ±0.91. Figure 5 shows the 68%, 95% and 99% CL allowed regions in the (S, T) plane 
for freely varying U (left) and the constraints found when fixing U = (right). For illustration 
also the SM prediction is shown. The Mh measurement reduces the allowed SM area from the 
grey sickle, defined by letting Mh float within the indicated range, to the narrow black strip. 



4 Conclusion 



Assuming the newly discovered particle at ~126 GeV to be the Standard Model (SM) Higgs boson, 
all fundamental parameters of the SM are known. It allows, for the first time, to overconstrain 
the SM at the electroweak scale and to evaluate its validity. The global fit to all the electroweak 
precision data and the measured Higgs mass results in a goodness-of-fit p-value of 0.07. Only a 
fraction of the contribution to the "incompatibility" stems from the Higgs mass, which agrees at 
the 1.3(j level with the fit prediction. The largest deviation between the best fit result and the 
data is introduced by the known tension between from LEP and An from SLD, predicting 
respectively a larger (by 2.5cr) and smaller (1.9<r) Higgs mass, and by R® for which an improved 
calculation increased the deviation from the measurement from previously 0.8<r to 2.4<r. 

The knowledge of the Higgs mass dramatically improves the SM predictions of several key observ- 
ables. The uncertainties in the predictions of the W mass, sin 2 d^ s , and the top mass decrease 
from 28 to 11 MeV, 2.3 • 10 -5 to 1.0 • 10 -5 , and from 6.2 to 2.5 GeV, respectively. The improved 
accuracy sets a benchmark for direct measurements, which has been reached (and surpassed) only 
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for the top mass. Theoretical uncertainties due to unknown higher order electroweak corrections 
contribute approximately half of the uncertainties in the M\y and sin 2 6^ ff predictions. 

The results reported in this letter depend on the validity of the assumption that the observed 
particle is indeed the SM Higgs boson. New physics may lead to deviations in the couplings, which 
also affect the global fit. The next round of experimental updates by ATLAS and CMS will lead to 
a more precise assessment of the new particle's properties and are expected with great excitement. 

We thank Louis Fayard whose questions have triggered the more detailed error analysis for the Mw and 
sin 2 6^ ff predictions provided in this version of the paper. 
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